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Table 2. W UMa-type Systems in Globular Clusters 



name P (days) V^a^ B-V V -I Amply M^*"" M^^ M^' Cl'^^ Cl^' Type 



NGC 104 





V214 


0.2737 


17.96 




0.82 


0.38 


4.59 




4.98 




1 






V221 


0.3135 


17.78 




0.79 


0.44 


4.41 




4.61 




1 




o , 


V225 


0.2346 


19.47 




1.04 


0.53 


6.10 




6.07 




1 




o ■ 


V227 


0.3788 


16.49 




0.52 


0.28 


3.12 




3.27 




1 


BS 


o ■ 


V238 


0.2506 


18.46 




0.77 


0.34 


5.09 




4.97 




1 




<N : 


V244 


0.3837 


16.16 




0.51 


0.22 


2.79 




3.20 




1 


BS 


a ■ 

■ 


V245 


0.2789 


15.49 




0.69 


0.38 


2.12 




4.47 




N 




V249 


0.3226 


17.33 




0.64 


0.33 


3.96 




4.01 




1 




1—5 ; 
o ■ 


NGC 288 
























VIO 


0.4388 


19.09 


0.61 




0.38 


4.40 


3.46 




2 








NGC 4372 






















152vl 


V4 


0.2418 


18.88 


0.83 


1.23 


0.38 


3.87 


4.19 


5.13 




N 




V5 


0.3403 


19.07 


0.83 


1.29 


0.42 


4.06 


3.53 


4.69 


2 


2 




V16 


0.3084 


17.84 


0.81 


1.02 


0.27 


2.83 


3.66 


3.90 


2 


N 


BS 


V22 


0.4150 


16.84 


0.63 


0.78 


0.14 


1.83 


2.54 


2.46 


2 


2 


BS 




V27 


0.3772 


17.48 


1.46 


1.76 


0.56 


2.47 


5.23 


6.20 


N 


N 




o : 


V31 


0.3777 


17.47 


1.12 


1.25 


0.45 


2.46 


4.20 


4.35 


N 


N 




o . 


V33 


0.2791 


19.09 


1.38 


1.54 


0.56 


4.08 


5.57 


5.98 


N 


N 




o ■ 


V35 


0.2901 


17.40 


1.19 


1.38 


0.69 


2.39 


4.92 


5.32 


N 


N 






Rup 106 

V7 


0.2755 


19.42 


0.98 




0.28 


2.17 


4.96 




N 






o ■ 


V17 


0.4135 


16.60 


0.84 




0.48 


-0.65 


3.76 




N 






NGC 5139 






















■ 


VIO 


0.3687 


17.31 




0.69 


0.34 


3.34 




3.57 




1 


BS 


iv: 


Vll 


0.3073 


17.20 




0.52 


0.30 


3.23 




3.31 




1 


BS 


V12 


0.2759 


17.79 




0.56 


0.31 


3.82 




3.66 




1 


BS, E 


><■ 


V13 


0.3055 


17.12 




0.44 


0.32 


3.15 




3.03 




1 


BS 


:ar 


V19 


0.3982 


16.22 




0.38 


0.29 


2.25 




2.30 




1 


BS 


V20 


0.3418 


16.60 




0.42 


0.16 


2.63 




2.74 




1 


BS 




V21 


0.2493 


17.86 




0.66 


0.19 


3.89 




4.22 




1 


BS 




V44 


0.2963 


17.26 




0.60 


0.23 


3.29 




3.67 




1 


BS 




V48 


0.3319 


16.37 




0.61 


0.11 


2.40 




3.49 




N 


BS? 




V49 


0.3663 


18.57 




0.92 


0.73 


4.60 




4.42 




1 


EB 




V54 


0.2838 


18.94 




1.09 


0.24 


4.97 




5.53 




2 






V56 


0.2812 


19.23 




1.12 


0.50 


5.26 




5.66 




1 






V57 


0.4178 


16.95 




1.03 


0.14 


2.98 




4.57 




N 






V58 


0.4158 


16.82 




0.64 


0.10 


2.85 




3.16 




1 


BS 




V61 


0.6309 


16.19 




0.52 


0.04 


2.22 




1.92 




1 


BS 




V64 


0.3851 


16.47 




0.48 


0.09 


2.50 




2.73 




1 


BS 




V65 


0.5122 


17.16 




0.61 


0.07 


3.19 




2.65 




2 


BS 




V67 


0.2506 


17.52 




1.70 


0.17 


3.55 




7.98 




N 





Table 2. (continued) 



name 


P (days) 


T r 

ymax 


B — V 


V — 1 


Amplv 


71 /fobs 
My 


My My 




CI 


Type 


NGC 5272 




















V238 


0.4900 


17.25 


0.58 




0.10 


2.13 


3.22 


N 






NGC 5466 




















NH19 0.3421 


18.42 


0.15 




0.28 


2.27 


2.64 


1 




BS 


NH30 0.2975 


19.23 


0.16 




0.32 


3.08 


2.94 


1 




BS 


NGC 


6121 




















V44 


0.2637 


17.69 


1.00 




0.23 


4.86 


4.62 


1 






V47 


0.2700 


16.85 


0.82 




0.29 


4.02 


4.03 


1 






V48 


0.2825 


16.42 


0.70 




0.22 


3.59 


3.58 


1 




BS 


V49 


0.2976 


16.95 


1.18 




1.02 


4.12 


4.93 


2 




EB 


V50 


0.2665 


17.17 


0.83 




0.44 


4.34 


4.09 


1 




EB 


V o ± 


U.OUOl 


17.03 


0.94 




0.42 


4.20 


4.17 


1 






V53 


3085 


15.75 


0.59 




0.25 


2.92 


3.08 


1 




BS 


V55 


0.3108 


16.70 


0.81 




0.49 


3.87 


3.73 


1 




EB 


NGC 


6362 




















V39 


0.3633 


17.84 




0.56 


0.23 


3.05 


3.31 




1 


BS, E] 


V43 


0.2850 


16.02 


0.73 




0.44 


1.23 


4.50 


N 






V45 


0.4090 


16.93 


0.76 




0.30 


2.14 


3.90 


N 






V52 


0.3990 


15.36 


0.64 




0.50 


0.57 


3.58 


N 






NGC 


6397 




















V7 


0.2716 


16.80 


0.61 




0.38 


4.44 


3.93 


2 






V8 


0.2710 


16.19 


0.50 




0.37 


3.83 


3.60 


1 




BS 


NGC 


6441 




















V48 


0.6674 


15.23 




1.01 


0.32 


— 1.39 


2.15 




N 


BS? 


V49 


1.0100 


16.52 




1.07 


0.36 


—0.10 


1.57 




N 


BS? 


V50 


0.4335 


17.85 




1.30 


0.55 


1.23 


4.04 




N 




NGC 


6752 




















V4 


0.2502 


18.49 


1.09 




0.55 


5.36 


5.96 


2 






V5 


0.2736 


16.89 


0.81 




0.52 


3.76 


4.94 


N 






V6 


0.3007 


19.21 


0.40 




0.45 


6.08 


3.52 


N 




BS? 


V8 


0.3150 


17.08 


0.29 




0.40 


3.95 


3.10 


2 




BS 


V9 


0.3630 


15.02 


0.61 




0.05 


1.89 


3.80 


N 






Vll 


0.2783 


15.91 


0.88 




0.30 


2.78 


5.12 


N 






V14 


0.3175 


18.78 


0.94 




0.13 


5.65 


5.05 


2 






NGC 


6838 




















VI 


0.3489 


18.29 




1.09 


0.50 


4.54 


4.55 




1 




V2 


0.3672 


17.81 




0.94 


0.65 


4.06 


3.90 




1 




V3 


0.3739 


19.07 




1.40 


0.35 


5.32 


5.54 




1 




V5 


0.4045 


17.94 




0.98 


0.45 


4.19 


3.86 




1 





*V4 in NGC 4372 relegated from Class-1 to Class-2 member because of the large uncertainties about the 
cluster distance modulus and reddening; see the text. 

Notes to Table 2. 

The data in the first six columns have been copied from discovery papers, as cited in Table 1. M^'" are 
the absolute magnitudes derived from distance moduli while My^ and My' are the corresponding absolute 
magnitudes derived from the respective calibrations. The columns labeled Cl^^ and Cl^' give the membership 

class, as explained in the text; N in this columns means that the system is not a cluster member. The last 
column contains the code EB for contact systems with EB-type hght curves and the code BS for Blue Stragglers. 



AMv(BV) 




Table 3. Expected effects of metallicity variations 



[Fe/H] 


A(B - V) 


A{V-I) 


Ai? 


Ai? 


AAfy 


AMy 








C = 5 


C= 10 


C = 5 


C== 10 


0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


-0.50 


-0.057 


-0.031 


-0.042 


-0.014 


0.280 


0.092 


-1.00 


-0.098 


-0.049 


-0.086 


-0.049 


0.572 


0.325 


-2.00 


-0.131 


-0.039 


-0.135 


-0.094 


0.900 


0.629 


-3.00 


-0.145 


-0.028 


-0.118 


-0.077 


0.785 


0.510 




0.65 



0,0 



-0,5 



-1,0 -1,5 
[Fe/H] 



-2,0 



(B - V)o 




(B - V)o 
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o 
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ABSTRACT 

A sample of 86 contact binary systems in 14 globular clusters with available color 
index data in [B — V) or in iy — I) has been analyzed. A large fraction of all systems 
(at least one third) are numerous foreground Galactic Disk projections over long lines 
of sight to the clusters. Since the selection of the cluster members has been based 
on the My (log P, coZor) calibrations, the matter of a metallicity-correction required 
particular attention with the result that such a correction is apparently not needed 
at the present level of accuracy. Analysis of the color-magnitude and period-color 
relations shows that globular cluster members are under-luminous relative to the 
Galactic Disk contact systems mainly because of the smaller sizes and, consequently, 
shorter orbital periods; the color-index effect of the diminished blanketing is less 
important, especially for {V — I). Among the Class-1 members (deviations in My 
smaller than 0.5 mag.), the most common are Blue Straggler (BS) systems. The 
variability amplitudes for the BS systems show a significantly different distribution 
from that for systems below the Turn Off Point (TOP): The BS systems in the sample 
have only small amplitudes while the distribution for the systems below the TOP is 
peculiar in containing only large amplitude systems. This difference may be linked to 
the relatively small number of the detected Main Sequence contact systems below the 
TOP as resulting from an observational selection effect due to the rapidly increasing 
measurement difhculties below the TOP. As a consequence, efforts at determining the 
frequency of occurrence of the contact systems below the TOP have been judged to be 
premature, but the frequency among the BS stars could be moderately well established 
at about 45 =b 10 BS stars per one contact BS binary. 

Subject headings: clusters: globular - binaries: eclipsing - blue stragglers 

1. Introduction 

Our thinking about close binary stars in globular clusters (GC's) has undergone a tremendous 
change in the last decade. Once such binaries were thought to be totally absent in the GC 
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environment, while now they appear to be of great importance to the dynamical evolution of 
the clusters. A large volume of research on the dynamical effects of the binary systems on the 
cluster evolution, including important and complex inter-relations with the evaporation and 



tidal-striping effects, has been summarized in several reviews starting with Hut et al. 1992, with 



updates in Sections 9.5 and 9.6 of [Meylan &: Heggie 1997] and in McMillan et al. 1998 . Solar-type 



contact binaries (also known as W UMa-type variable stars) contain the least amounts of angular 
momentum that binary systems made of Main Sequence components can have. They represent 
last stages of the angular momentum loss (AML) evolution of primordial binaries or are one of 
the products of the dynamical inter-cluster interactions. Together with the RR Lyr-type pulsating 
stars on the horizontal branch and SX Phe-type pulsating stars among Blue Stragglers, the 
W UMa-type binaries are the most common type of variables in the GC's. 

The number of W UMa-type binaries in GC's has expanded in the recent four years from 24 
( Mateo 1996a| ) to the current number of 86. This paper attempts to integrate and analyze the 



available data for these variables in order to establish most essential results and to guide further 
research. The analysis is similar in its goals to the study of the Galactic Disk systems which 
compared the W UMa binaries in the galactic field - as seen in the direction of the galactic Bulge 
in the OGLE survey - with those in several old open clusters ( Rucinski 1998b| = R98). Searches 



for contact binaries in GC's are much more difficult than in open clusters. They require - on top 
of a generous allocation of observing time permitting variability detection and monitoring over 
several nights ~ at least moderate-size telescopes located in sites with excellent seeing. Cores of 
some clusters remain too compact for photometry of individual stars even in perfect seeing; such 
clusters must be observed with the Hubble Space Telescope for which continuous monitoring of 
variability over long periods of time is difficult to arrange. It is expected, however, that new 



methods of analysis of difference images, such as developed by Alard fc Lupton 1998 and Alard 



1999, will result in substantial progress in detection and analysis of variable stars in very dense 
fields close to the cores; one of the first applications to GC's (RR Lyr stars in the core of M5) is 
by lOlech et al. 199E . 



Because of the action of two possibly mutually reinforcing types of processes, of dynamical 
interactions and of the magnetic braking, rather than one - as in the Galactic Disk and in open 
clusters, where only the latter mechanism would be sufficiently effective - one may expect high 
frequency of occurrence of the contact systems in globular clusters, possibly even higher than in 
the Disk. It is, in fact, surprisingly high in the Disk: As shown in R98, the Disk systems show an 
increase in the frequency of occurrence among F-K dwarfs over time, in the accessible range of 0.7 
- 7 Gyr, reaching the spatial frequency as high as about one such binary (counted as one object) 
per 80 - 100 single stars in the Galactic Disk. Contact systems with spectral types earlier than 
about middle A-type and orbital periods longer than 1.3 - 1.5 days are less common with currently 



un-measurable frequency of occurrence, as shown in Rucinski 1998a; this long-period cutoff may 



be a function of the parent population because the contact binaries with periods up to 2.5 - 3 



days are known to exists in LMC ( Rucinski 1999 ). The high frequency of contact binaries with 
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a gradual increase with age in old open clusters is consistent with the prolonged magnetic-wind 
braking AML process acting over a time scale of some 1-5 Gyr (R98) and producing relatively 
long-lived contact systems. 

Contrary to expectations based on the above reasoning, the preliminary indications from 
individual studies that reported discoveries of contact binaries in GC's (and are cited in this 
paper) do not confirm the high frequency of occurrence in these clusters. As several authors of 
such papers already remarked - but bearing in mind the tremendous technical difficulties - the 
frequency appears to be relatively low, at the level of a small fraction of a percent and some 
clusters do not seem to have W UMa-type binaries at all. An attempt to assess this matter is 
presented in this paper. It is argued that, at this moment, we cannot really say much about the 
frequency of occurrence of the contact systems on the Main Sequence, below the Turn-Off Point; 
however, the frequency among the Blue Stragglers appears to be high, some 2 to 3 times higher 
that among the stars of the Galactic Disk. 

The current paper consists of a description of the sample of clusters with contact binaries 
in Section |2|, then of the sample itself in Section ^. The observed metallicity effects affecting 
the absolute-magnitude calibration (which is used to select the members) and affecting the 
observed properties of the systems are described in Sections ^ and ^. The color-magnitude and 
the period-color relations are discussed in Sections |^ and 0. Details concerning systems with 
EB-type light curves. Blue Straggler contact systems and the frequency of occurrence are given in 



Sections q - 10. Conclusions are stated in Section 11 



2. The cluster sample 

The sample of GC's surveyed deep enough to include Main Sequence stars, below the Turn-Off 
Point (TOP), currently consists of 14 clusters. The sample is quasi random in the sense that 
several authors contributed the data using their own preferences, but that practically all data 
have come from ground-based telescopes. Dr. Kaluzny and his collaborators, who contributed 
most of the results, selected primarily the nearest clusters, with moderately developed cores 
(permitting photometry close to the centers) and avoiding those with small galactic latitudes, 
within |6| < 10°. The essential parameters characterizing the clusters are given in Table |l|, with 
clusters arranged according to the NGC number. To insure uniformity of these parameters, they 
have been taken from the database of parris 1996 , version June 22, 1999, which is available at: 



http: / /physun. physics. mcmaster.ca/Globular. html . In Table || we give the following parameters: 



the galactic coordinates I, b in degrees, the galacto-centric distance Rgc in kpc, the reddening 
Eb-Vj the observed distance modulus (m — M)v, the metallicity parameter [Fe/H] and the 
concentration parameter c (for collapsed cores, c = 2.5). The galacto-centric distances span 
a wide range 3.5 < Rgc < 18.5 kpc while the metallicities occur within the representative 
range —2.22 < [Fe/H] < —0.73. The sample is dominated by clusters with moderate and low 
concentration; clusters with c 2.5 are under-represented. Eight clusters have galactic latitude 
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|6| < 20° indicating a possibility of large interstellar extinction and heavy Milky Way stars 
contamination; for three clusters, NGC 4372, NGC 6121 and NGC 6441, the reddening is large, 
Eb-v > 0.3. 

The cluster sample requires some comments: 

1. While Table ^ lists 14 clusters, the entry for NGC 5904 (M5) actually reports a null result 
because all the systems suggested as contact binaries by |Yan &: Reid 1996 turned out to 



be spurious detections, as has been shown by Kaluzny et al. 1999| . The cluster has been 



monitored for variability of its stars by several investigators (see references to Table ||). 

2. NGC 4372 is a very important cluster because of its low metallicity and a large number 
of systems discovered in its direction. Unfortunately, because of its low galactic altitude 



of b = —10°, it is has a large and patchy reddening. The discovery paper by Kaluzny fc 



Krzeminski 1993 gives, in its Table 3, the observed data in (B — V) and {V — I), but the 
patchiness-corrected data are given only for the former color index. The corrections have been 
calculated for (V-I) assuming the relation AEy-i = 1.24:AEb-v = 1-24 [{B-V)c-{B-V)]. 
The same slope is used throughout this paper in deriving the values of Ey-i from Eb-v- On 
top of the patchy reddening, the values of the mean reddening and of the distance modulus 
are very uncertain for this cluster. While Kaluzny & Krzeminski assume Eb-v = 0.48 
and (m — M)v = 14.8, the Harris database (Table || here) quotes Eb-v = 0.39 and 
(m — M)v = 15.01. For consistency, the Harris set has been used here, but ~ as we discuss in 
the next Section - this choice makes a very important difference in ascertaining membership 
of systems detected in this cluster. 

The color-magnitude of the core of NGC 6752 has been studied with the Hubble Space 
Telescope ( pR,ubenstein &: Bailyn 1997] ). It indicates a relatively high frequency of binary 
stars at a level of 15 — 38 percent in the inner core, but below 16 percent beyond the core. 
The relative frequency must be therefore a strong function of the radial distance from the 
cluster center. Unfortunately, searches for variable stars in the clusters analyzed in this 
paper are not uniform in this respect: Some clusters were observed with the inclusion of 
the cores, some were observed only at some distance from the center where crowding was 
assumed to be tolerable. This casts a large uncertainty on any considerations involving 
numbers of contact systems and their frequency of occurrence relative to other stars. 



3. Cluster members 

All contact binary systems discovered in globular clusters are listed in Table |2|. The 
Fourier-analysis of the light curves was not used to verify the W UMa-type or the EW shape of 
the light-curves (as in the selection of the OGLE sample in |Rucinski 1997a ) because of the partial 
coverage of some of the light curves which could produce incorrect values of the Fourier coefficients; 
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the general appearance of the light curve and the original classification by the discoverers were 
the only criteria used here. Contact systems with unequally deep eclipses (EB-type light curves) 
have been retained; they are marked as such in Table |^. Systems having light curves suggesting 
detached components (EA type) are not considered here. 

The empty entries in Table ^ are due to the fact that most of the photometric searches have 
been done either in B and V band-passes or V and Ic band-passes (the subscript indicating that 
the / band is of the Cousins system will not be used from now on). Thus, the available data split 
in two sets, forcing us to discuss all relations and all properties in parallel in V and (B — V) and 
in V and {V — I). We will call these the BV-set and the Vl-set. Only one cluster was observed 
in all three bands, NGC 4372 ( [Kaluzny fc Krzeminski 1993| ) permitting consideration of two 
photometric indices, {B — V) and {V — I), for 8 contact systems in this cluster. As luck would 
have it, this is the cluster which is the most heavily reddened with a highly patchy and uncertain 
extinction. 

Table || lists the variables with the names and designations as assigned in the discovery 
papers. For each system, we give the orbital period in days, then the maximum brightness 
V, {B — V) and/or {V — I) and the total variability amplitude Ay. The last columns give 
our results on the membership to the parent clusters and the variability /membership type (see 
below). Assignment of the membership has been performed by comparison of "observed" absolute 
magnitudes, derived from the distance modulus: My"^ = V — [m — M)y, with those derived from 
the calibrations, M^^ = My (log P, 5 - V) or = My (log P, V - I), as described in more 

detail below. A simple criterion for membership was used calling systems with the deviations in 
AMv = M0>' - M^'^' smaller than 0.5 mag the Class 1 members and those with 0.5 < \AMy\ < 1.0 
the Class 2 members (here My"^ takes the meaning of My^ or My^ depending which one has 
been available). Deviations larger than one magnitude were assumed to signify that the binary is 
not a member but a foreground or background projection. 

The calibrations used in this paper have been: 

M^^ = -4.44 log P + 3.02 {B - V)o + 0.12 (1) 
M^^ = -4.43 log P + 3.63 {V - I)o - 0.31 (2) 



The My^ calibration is based on the Hipparcos data ( Rucinski fc Duerbeck 1997| ) while the one 



for My^ was developed and served well for the analysis of all the OGLE data ([Rucinski 1997a 



Rucinski 1997b , Rucinski 1998a , Rucinski 1998b| ). No allowance for lowered metallicity of GC's in 



the calibration has been made. We discuss this important issue in the next Section 

The deviations AMy are shown in graphical form in Figure [l| versus the orbital period and 
color index, for both color-index sets. While many systems do fall within the band of AMy 
close to zero, a large fraction of systems are actually foreground projections onto the fields of the 
observed clusters. This fully agrees with the high frequency of the W UMa-type binaries in the 
Galactic Disk. The data based on the {V — I) index show a somewhat better consistency, with 
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smaller scatter in AMy and with a better definition of the group of the Class-1 members. This 
may be due to a weaker dependence of this index on the interstellar reddening and/or on the 
metallicity, but may possibly be simply due to lesser photometric difficulties of observing red stars 
in the V and / bands than in the B and V bands. As pointed out to the author by Dr. Kaluzny 
(private communication), selection of the best photometric bands is not an easy matter so that 
{V — I) may not always be preferable to {B — V): The upper MS is better observable in B than 
in / because the numerous red dwarfs from the lower MS produce a strong background in /; the I 
band is also inconvenient for red stars because Asymptotic Giant Branch stars are usually strongly 
over-exposed and prevent photometry of other stars close to cluster centers. 

The two penultimate columns of Table |l] give the number of contact systems in the GC's. ridet 
is the total number of systems in a cluster, while uqi and nc2 are the numbers of Class-1 and 
Class-2 members, respectively, nss is the number of Blue Stragglers (of both classes). Nbs are 
very approximate estimates of Blue Stragglers which were monitored for variability in the GC's. 



We will discuss these data later, in Sections and 10. The approximate locations of the color 
indices at the Turn-Off Point are marked by vertical broken lines in Figure |l|. Note that practically 
all systems to the blue of these lines (i.e. the Blue Stragglers) are members of the clusters. 

The cluster-member selection process described above assumes that the calibration formulae, 
Eq. H and Eq. ^ are applicable to contact systems in globular clusters. This is verified by assuming 
that most typical systems which have been detected in a given direction are actually genuine 
members of the clusters. In effect, we require that the calibrations reproduce the modal (most 
probable) values of the deviations AMy . The histograms of the deviations are shown in Figure ^ 
Disregarding a complication of the metallicity-dependence of My^'^ which will be discussed in the 
next section, we can observe the following facts: 



1. There are 35 Class-1 systems in the sample, that is roughly 1/2 of the total. 

2. The number of the cluster members would increase by four if we add Class-2 systems which 
have blue color indices and short periods characteristic for Population II Blue Straggler 
systems (Section |9|). 

3. We have no good argument to claim that any of the eight remaining Class-2 systems is a 
cluster member. 

4. There are 21 definite non-members in the sample, that is about 1/3 of the total number. All 
but one are foreground Disk systems. 

5. The only obvious background system is V6 in NGC 6752. The AMy deviation of about 2.5 
mag. suggests that the system is some 3 times further away than the cluster, at a distant 
periphery of the Galaxy at the distance of some 15 kpc. 



6. 



As discussed in Kaluzny et al. 1998a, there exists an ambiguity with the orbital period for 
the only contact system detected in M3, V238. It appears that neither of the acceptable 



-7- 



periods places the system in the cluster. The two possibilities are joined by a dotted line in 
Figure |^. Dr. Kaluzny (private communication) suspects - from analysis of the individual 
CCD images - that blending with a red giant is the cause of photometric difficulties with 
V238. 



4. Metallicity effects in the My calibrations 

The cluster-selection process requires a clear conceptual separation of the effects of lowered 
metallicity on the simplified My calibrations given by Equations |l| and ^ (they are called here 
M{f' or specified as M^^ or M^) which we use for selecting the cluster members from the 
effects genuinely influencing the observed binary properties, such as the absolute magnitudes 
My^ (derived from cluster moduli), the de-reddened color indices or the orbital periods. Here we 
concentrate on the effects affecting My^ and solely from the observational point of view, without 
going into the details of how genuine properties are affected; this will be discussed in the next 
Section |5|. We note that once we select a sample of the cluster members, we will discuss only the 
directly observed quantities; the values of My^ will not be used from this point at all. 



In Rucinski 1995 arguments have been presented that the My^ = My (log P, color) calibrations 



require small, but significant corrections for metallicity variations. The corrections would reflect 
the blue colors of Population II stars which have low atmospheric blanketing due to weak spectral 
lines of metals. It was argued that while the period term in Equations || and ^ would account for 
any differences in the system size, the color-index term would require a correction as a proxy of 
the effective temperature. The data and the available calibrations were very preliminary at that 
time; the current material is much richer so that the need for the corrections can be re-evaluated. 
This matter is very closely linked to the selection of the cluster members, so that we must explain 
the details of this process. 

Figure |3| shows the same deviations AMy as in Figure but this time plotted versus [Fe / H] 
for the clusters. The slanting lines give the metallicity corrections SMy^ = — 0.3 [Fe/H] and 
5M^^ = -0.12[Fe/H], as suggested in [Rucinski 1995| . To visualize how the corrections would 



impact the membership selection process, the AMy deviations - with the metallicity corrections 
applied - are shown as broken-line histograms in Figure 

Figures |2| and |3| contain the only available information on the need of the corrections and how 
their use would influence the selection of the cluster members. The situation is relatively simple for 
the VI set: Because of the well-known low sensitivity of this color index to metallicity variations, 
the current data do not forcefully suggest or reject the need for the metallicity corrections, 
although the case of no corrections seems to be slightly preferable; for simplicity, we assume that 
they are not needed. The data for the BV-set show a large scatter in the band of the expected 
members and its is difficult to decide if the corrections are really needed. The crucial clusters 
for resolution of the problem for the BV set are the two clusters with the lowest metallicities 
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[Fe/H] < -2, NGC 4372 and NGC 5466. 

As remarked before, NGC 4372 is seen at low galactic latitude through a patchy interstellar 
extinction. However, the cluster was observed in BVI so that two sets of color indices are available 
and a consistency check is in principle available. Assuming the data as in Harris' database 
(as given in Table |l|), the My^ calibration indicates that none of the contact systems can he 
classified as Class- 1 member while the My^ calibration suggests that one system, V4, is a Class 1 



member. However, if we follow the assumptions of [Kaluzny fc Krzeminski 1993| , Eb-v = 0-48 
and {m — M)v = 14.8, then large shifts in AMy by 0.49 and 0.61 mag. for both color-index sets 
occur: While M^^" become fainter by +0.21 mag., the predicted My^ and My^ become brighter 
by —0.28 and —0.40 mag. As the result, all data points for this cluster slide down in Figure ^ 
by amounts shown there by arrows. While for the Harris data only the system V4 would be a 
Class- 1 member in the BV set (without a confirming evidence from the VI set), now the systems 
V4, V16 and V22 would be Class-1 members for the BV set and V5, V16 and V22 would be 
Class-1 members for the VI set. One gains then in consistency between both sets for V16 and 
V22, but the matter of membership remains unclear for V4 and V5. It is obvious that only three 
or four systems among the eight might be members of NGC 4372, but we cannot be absolutely 
sure which ones. In this situation, we have taken a conservative approach and conclude that the 
data for NGC 4372 are too uncertain to be sure of the membership of the systems in this cluster; 
this cluster cannot tell us much about the need for a metallicity correction in the calibration. We 
note that judging by its period-color combination, the system V22 is a genuine Blue Straggler and 
irrespectively how big is its deviation AMy, it is almost certainly a member of NGC 4372 (see 
Section |9|) . 

NGC 5466, with even more extreme metalhcity than that of NGC 4372 of [Fe/H] = -2.22, 
follows the Hipparcos calibration for the Disk stars ( Rucinski fc: Duerbeck 1997] ), without any 



metallicity correction, very well. The two contact binaries must be genuine members simply 
because at the galactic latitude of 6 = +74° chances of having Disk population stars within the 
cluster are practically equal to zero. Also, both stars are exceptionally blue, belonging to the Blue 
Straggler group of contact binaries (see Section P); equally blue systems very rarely occur in the 
Galactic Disk. 

Guided mostly by the case of NGC 5466 and by the resulting simplicity of the assumption, 
we conclude that the My calibrations established for [Fe/H] = apparently work well for low 
values of metallicity and that the [Fe/ff] -corrections in the expressions for M^"' are apparently 
not needed. As we will see in the next sections, the color indices of the low-metallicity systems 
are definitely much bluer than those for Disk systems, so that - without any corrections ~ 
one would expect values of My^ indicating artificially higher luminosities. Since this is not 
observed, some other metallicity-dependent factor must provide a compensating effect through 
the period-term (this may be for example the un-accounted for influence of the mass-ratio). 



Concerning the previous, apparently erroneous result in Rucinski 1995, we remark here that the 



correction of —0.3 x [Fe/H] for the [B — y)-based calibration was suggested for an old version 
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of the M^^ calibrat ion which has been superseded by the much better Hipparcos caHbration. 
The corresponding correction for the {V — I)-caHbration of —0.12 x [Fe/H] was suggested for 
consistency with the one for {B — V), but it was always recognized that it was smaher than the 
measurement and definition uncertainty of the cahbration itself. 



5. Expected effects of low metallicity 

In summarizing the main effects of lowered metallicity, we discuss two effects which manifest 
themselves very differently: of the blue atmospheres (for the same effective temperature) and of 
smaller stellar sizes. 

The decreased metallicity influences the atmospheric structure in that metal lines produce less 



blanketing so that stars become bluer. Only this effect was discussed in [Rucinski 1995| , together 
with its possible influence onto the My calibrations. For consistency with the previous results, we 
use the same relations between the color- index changes and [Fe/H] as evaluated from the models 



of Buser &: Kurucz 1992| . The expected changes in the color indices {B — V) and {V — I) are given 



for a MS star atmosphere at Tg// = 5000 K in Table |3|. Obviously, several qualifications may be 
in order here: Single-star model-atmosphere results may be in-applicable to magnetically active 
contact binaries and the color index data have been calculated for spherical stars whereas strong 
and variable limb darkening effects are always present in contact binaries. 

The low-met allicity systems are expected to have smaller dimensions and - for the same 
contact geometry as for Galactic Disk systems - should have shorter periods than Population I 
systems. In the currently sole study on properties of Population II contact systems binaries by 



Webbink 1979| , the stress was on the influence of the prolonged angular-momentum loss for such 
old objects. The author, however, remarked about the smaller sizes of such systems relative to 
Population I systems, but did not discuss this point any further. Some guidance on the expected 
effect can be found in the mass-radius relation and its dependence on the stellar metallicity. For a 
fixed mass, Kepler's law enforces the proportionality logP cx 3/2 log A. If the geometry of contact 
is the same irrespective of the chemical composition so that the relative sizes of components are 
independent of the metallicity, rj = Ri/A, then smaller stellar sizes R should lead to smaller 
A and to shorter orbital periods. Again, a qualification may be in order: The period changes 
may depend on metallicity in a much more complex way than just through simple radius scaling 
because the internal structure of Population II contact binaries does not have to be identical to 
that of Population I systems. 

The mass-radius relation for low- met allicity stars is currently a subject of very lively 
discussions, stimulated by the distance determinations for subdwarfs, as provided by the Hipparcos 
mission. One of the most recent theoretical investigations of stellar models with varying 
abundances is by pastellani et al. 1999| . The absolute magnitude calibration for stars with 
Tgff = 5000 K (Eq. (1) in this paper) can be re-written into a metallicity-radius dependence: 
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Alogi? = 0.158 [Fe/H] + 0.035 [Fe/Hf - 0.485Ay; with the solar abundance assumed to be 
Zq = 0.0169 and [Fe/H] = \ogZ — log Zq. The size of the correction for hehum abundance 
changes, AY, which accompanies the metalhcity changes as the stellar population ages, is a 
difficult matter. The authors discuss the large uncertainty in the ratio C = /^Y/ /\Z which is 
currently very poorly known; some observational results suggest C = 3 it 2 while the authors 
consider C ~ 5 — 6. The metallicity-radius relation is shown in Figure ^ for three values of C. 

With the expected values of A log P oc 3/2 x A log R one can compute the expected variations 
of the absolute magnitude. We have a choice here: We can simply assume that the luminosity 
will scale with the square of the radius or we can use the log P terms in the Equations |l] and ^. 
The former choice would suggest: AMy = —5 x Alogi?, while the latter choice gives a steeper 
dependence: AMy = —4.44 x 3/2 x Alog-R (the period-term coefficients are basically identical 
for both color- index sets). In what follows, we will use the second expressions for consistency with 
the adopted expressions for My\ observing that the steepness of the period term (which probably 
hides many unaccounted period-dependent effects) may actually explain the unexpected absence 
of the [Fe/H] term in the equations giving MyK 



6. Color-magnitude diagram 

The color-magnitude diagrams for the sample of the GC members are shown in Figure ^. The 
data plotted are the observed My"^ and the de-reddened color indices so that uncertainties with 
the metallicity corrections in My^ do not enter directly into this figure (but only through the 
selection of cluster members, via the deviations AMy). The Class-1 members, which we consider 
genuine members of the clusters, are marked by filled circles and the range of metallicity {[Fe/H] 
smaller or larger than —1.5) is shown by the size of the symbol. The figure also contains the 
theoretical isochrones computed by the Padova group (iBertelh et al. 1994]) for 14 Gyr and for 



three values of metallicity [Fe/H] = —1.66, —1.26 and —0.66. In addition, the data for members 
of old open cluster members are shown by small x-symbols, following the results in R98. 

The most striking feature of the color-magnitude diagrams for Population II contact binaries 
is the shift of the contact-binary sequence by about one magnitude below that for Galactic Disk 
systems. A similar shift is well known for single subdwarfs, but it is seen here for contact binaries 
for the first time. The sequence is relatively well defined and extends uniformly on both sides of 
the Turn-Off Point (TOP), similarly as in old open clusters such as Cr 261 (R98). 

A group of Blue Stragglers (BS) in the low-metallicity clusters, to the blue of the TOP is the 
most conspicuous group of the contact binaries in the color-magnitude diagrams in Figure |5[ A 
few Class-2 systems also have very blue color indices of the BS stars; since such systems practically 
do not occur among Disk contact systems we think that they are not foreground projections, 
but genuine cluster members with the AMy deviations larger than 0.5 mag. These systems are 
in the BV set: V22 in NGC 4372 and V8 in NGC 6752; in the VI set: V22 in NGC 4372 and 
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ves in NGC 5139. Thus, V22 in the controversial cluster NGC 4372, with its Blue Stragglier 
characteristics, is intrinsically the brightest system in the curr ent GC sample with M^'" = 1.83, 
(B - V)o = 0.24 and (V - I)o = 0.30. 



7. Period-color diagram 

The period-color diagram is, to some degree, a more natural and precise way to display 
properties of contact binaries than the color-magnitude diagram. This is because one quantity, the 
orbital period, is known basically without error, at least when compared with photometric errors. 
Only one photometric quantity - the color index - then enters into the picture. The color index 
can still be affected by several factors, in addition to the measurement errors, the most severe 
being the uncertainty in the reddening correction, but the period-color relation is not affected by 
errors in the distance modulus. 

Figure |^ shows the two period-color diagrams for the GC sample, with the same symbols as 
in Figure |5[ The reader is suggested to view both figures simultaneously and to note the common 
features. The figure contains the Short-Period Blue-Envelopes (SPBE) for Disk systems, shown 
as continuous lines for both color-index sets ( [Rucinski 1998b| , Rucinski 1997a| ). Their shapes are 



given by: {B - V)spbe = 0.04P-2-25 and {V - I)spbe = 0.053 P-^-i (the orbital period P is 
in days). While the numerical values in these definitions do not have any physical meaning, the 
curves are important because they delineate location of the least-evolved contact systems in the 
Disk Population sample. As many theoretical investigations indicated, contact systems would 
normally evolve away from the main Sequence, toward larger stellar dimensions, longer orbital 
periods and cooler atmospheres. The interstellar reddening also increases the color index. Thus, 
the SPBE has a meaning of a Zero- Age Main Sequence for the Disk population systems. 

The above considerations should also apply to Population II contact systems except that, as 
pointed out by Webbink, the angular-momentum loss due to gravitation radiation emission may 
win over the long time scales involved and eventually lead to shortening of the orbital periods. As 
we see in Figure |^ the GC systems indeed have periods shorter than those for Disk systems, but 
this is expected irrespective of whatever mechanism of forming them is involved: As discussed in 
Section |5|, to be a contact system, a low-metallicity binary must have a short period because its 
components are small. Judging by the size and direction of the metallicity corrections (shown by 
arrows), the relatively larger effect is observed in AlogP than in the color-index shifts A(i? — V) 
and A(y — /). Although the two changes can partly compensate each other in their control of 
the absolute magnitude, the compensation is not exact because the arrows in Figure ^ are not 
parallel to the lines of constant My ■ This was visible in the color- magnitude diagrams in Figure |^, 
where the GC systems were obviously fainter than the Disk systems. It is striking how different 
Population II contact systems are from the Galactic Disk systems, but also how similar they are 
within their group. Apparently, a difference in metallicity from the solar [Fe/H] ~ to —0.7 or so 
produces a relatively larger change in their properties than a further change to —2.2 observed for 
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the most extreme-metallicity systems. Since most changes are due to the change in the component 
dimensions, this would argue for a relatively small value of the currently poorly-known coefficient 
C = AY/AZ (see Figure | in Section |). 



8. EB-type systems 

The sample of binary stars considered in this paper contains six systems with EB-type light 
curves. Such light curves are characterized by unequally deep eclipses, but with strong variations 
between minima suggesting possibility of a physical contact. The OGLE sample ( Rucinski 1997b| ) 



contained only 2 systems of this type among 98 systems in a volume-limited sample of contact 
binaries selected using a Fourier light-curve shape filter. In this paper, we do not use this filter 
because it is sensitive to the phase coverage and tends to be too discriminatory. Thus, we assume 
that the six systems are related to contact systems. In fact, these may be various forms of 
semi-detached binaries in the pre-contact or broken-contact stages, with either the more massive, 
hotter or less-massive, cooler components filling their Roche lobes ( [Eggleton 1996| ). The latter 
cases appear to be less common, but a good case of a very short period Algol has been recently 



identified in W Crv (Rucinski & Lu 199 



Three among the six systems appear in the BV set and three in the VI set. Inspection of 
Figure || and Table |2| shows that all but one are Class-1 cluster members; the one with a slightly 
larger AMy is a Class-2 member. Thus, they all follow the absolute-magnitude calibrations for 
normal contact binaries and seem to be genuine cluster members. This relatively high frequency 
of occurrence among contact systems, of 6 among 35 (or 39 if 4 Class-2 Blue Stragglers are added) 
is unexplained and interesting. We note that two among the six systems are Blue Stragglers. 



9. Amplitude distribution and the Blue Stragglers 

It has been noted in Section ^ that the contact system sequence continues without any obvious 
changes in the period, color or absolute magnitude properties across the Turn-Off Point, into the 
Blue Stragglers domain. Now we will address the only property which is available to characterize 
the light curves, the variability amplitude and its distribution. 

The Blue Stragglers (BS) are an important group of stars in old stellar clusters. It is now 
recognized that they must form through binary evolution processes, although it appears that there 
are actually many such processes and it is not easy to find out which one occurs most commonly 
( Leonard 1996| , Mateo 1996b| ). The BS formation and evolution is such a large and active area 



that special meetings have been devoted to it ( Saffer 1993| ) and very active research continues. 



Contact BS are relatively easy to identify for several reasons: (1) They are bluer than typical 
galactic Disk systems which start appearing at {B — V)o < 0.3 (R98), (2) Their photometry is 
relatively less difficult than for the Main Sequence stars because they are photometrically well 
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above the level where - usually formidable - crowding problems for the MS stars set in. The 
question is: Are they in any other way different from the "normal" Main-Sequence contact systems 
in the GC's? One such property can be amplitudes of light variations. Although the amplitude 
statistics involves a convolution of the distribution of the mass-ratios with the distribution of 



orbital inclinations, as was discussed in Rucinski 1997b, lack of large amplitudes must mean that 



the large {q ^ 1, q = M2/M1 < 1) mass-ratios do not occur: When components differ in sizes, 
only small amplitudes are possible. 

Figure shows the observed amplitudes for the GC members plotted in relation to the 
intrinsic color indices. For the color index values, we assume that systems with (B — V)o < 0.4 
or {V — 1)0 < 0.55 are Blue Stragglers. The amplitudes do show a change at the Turn-Off 
Point, but this change is well defined only for the VI set: The large amplitudes are observed 
only for the systems to the red of the TOP, that is for the genuine MS systems. The two-sided 
Kolmogorov-Smirnov tests comparing distribution on both sides of the TOP's, limited to Class-1 
systems, indicates that the difference in the distributions is not significant for the BV set, but the 
probability of a random chance producing the observed difference for the VI set is only 1.1 x 10^^. 
The significance changes only slightly when Class-2 members are added with the probabilities 
of a random result of 0.05 for the BV set and 0.71 x 10~^ for the VI set. The change in the 
observed amplitude distribution at the TOP is therefore very highly significant for the VI set, but 
insignificant for the BV set. At this moment, we have no explanation for this difference between 
the color-index sets except that - as pointed several times in this paper - it may be related to the 
photometric difficulties in the i?-band for red stars and to the larger and more uncertain reddening 
corrections in the {B — V) color index. 

On the basis of the above numbers one would be tempted to claim that, indeed, the BS's 
have smaller amplitudes than the Main-Sequence contact systems. However, a caution is in order 
here: It is very important that we do not see small amplitudes among the Main Sequence systems 
below the TOP, only the large ones. They are obviously missing there because for random orbital 
inclinations systems showing small amplitudes should be always more common than systems 
showing large amplitudes. The histograms in the upper parts of the panels of Figure |^ contain 
the expected amplitude distributions calculated in [Rucinski 1997b for two assumptions of the 



mass-ratio distribution, a flat Q{q) and Q{q) = {1 — q). The data for the Disk systems in the 
OGLE sample suggested that the distribution may be actually even more strongly decaying with 
5 — > 1 than the one oc {1 — q), because large amplitudes are exceedingly rare in an unbiased 
sample (this is very much unlike the sky- field sample). In the case of the present GC sample, a 
comparison of the theoretical distributions with the observed ones suggests that the observed data 
may be severely modified by strong detection selection effects for systems below the TOP. While 
almost all systems in the BS group appear to be detected (there is only a small depression in the 
distribution for Amply < 0.2), the MS systems appear to be entirely missed for Amply < 0.3 due 
to the difficulties with accurate photometry in conditions where measurement errors and crowding 
problems rapidly increase with the apparent magnitude. 
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10. Frequency of W UMa-type systems in globular clusters 

The discovery selection effect described in the previous section may lead to an under-estimate 
of the number of the contact Main Sequence (below TOP) systems by a factor of the order of 
about 5 to 10. The previous results for the sky-field sample should be recalled here as a sobering 
experience: For several years the frequency of contact binaries of one per one thousand MS stars 
has been considered a well-established, "textbook- level" fact, in spite of the warnings that the sky 



sample of known contact binaries contained only large-amplitude systems (see Figure 2 in Kaluzny 



fc Rucinski 199^ ). Only the systematic characteristics of the OGLE sample (R98) have shown 



that the apparent]^ frequency is about ten times higher, reaching about 100 - 130 normal F-K 
dwarfs per one contact system. A very similar situation re-emerges here, except that this time we 
can directly suspect - from the difference in the amplitude distributions on both sides of the TOP 
- that the discovery selection effects are more severe below the TOP. We cannot correct for these 
selection effects because they must be different for each cluster and are very difficult to quantify. 

In this situation, it has been felt prudent to abandon attempts of determining the frequency 
below the TOP and concentrate on the frequency data for the BS systems. This is in turn 
complicated by lack of information in discovery papers on the number of Blue Stragglers which 
were monitored for variability. A simple, but potentially risky assumption has been made at 
that point that the numbers of the BS's shown in diagrams in the discovery papers are equal 
to the numbers of systems which were actually monitored. Frequently stars without measurable 
color indices (and not shown in color-magnitude diagrams) are monitored for variability, but one 
may hope that this would not happen in the BS region. Thus, approximate numbers of the Blue 
Stragglers, Nbs, have been estimated by simply counting of data points on the color-magnitude 
diagrams. It is stressed that these are very approximate estimates; for example, for uj Cen the 
available estimate for the fields D-F ( [Kaluzny et al. 1997b ) was multiplied by two to obtain the 



total number in all observed fields. For 47 Tuc ( Kaluzny et al. 1998b| ) unpublished data have 



been supplied by Dr. Kaluzny. The estimates are given in Table |^ and should be compared with 
the number of the contact BS systems, ubs- In doing so we immediately see the problem of 
numerically very small values for ubs which are obviously subject to relatively large Poissonian 
fluctuations. Thus, we can see that we are still very far from being able to correlate the frequency 
of contact BS systems with cluster properties such as the metallicity index, [Fe/H], or the 
concentration parameter, c. In this situation another disputable step was made in assuming that 
the average frequency is the same for all GC's. This can be derived by simply summing Nbs and 
Ubs for all clusters and taking their ratio. The result is that 20 contact binaries are observed 
among about 900 Blue Stragglers. Thus the average inverse frequency is Jbs = 45 it 10 normal BS 
stars per one contact system. This frequency is the apparent one, that is it applies to systems 



^We distinguish the apparent frequency, which is uncorrected for missed low-inchnation systems, from the spatial 



frequency which is about 1.5 to 2 times higher. The correction factor depends on the mass-ratio distribution ( Rucinski 



19971:) 
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which can be discovered, without any corrections for systems missed because of the low orbital 
inclinations. Since we do not know the mass-ratio distribution for the contact BS systems, we 
cannot correct for these missed systems to evaluate the true spatial frequency. 

The inverse apparent frequency of fss = 45 is significantly different from the one observed 
for the Disk stars, which is approximately foisk — 100 — 130. Thus, the Blue Straggler population 
of the globular clusters contains some 2 — 3 times more contact binaries than the Old Disk stars. 
We note also that the inverse frequency at the level of 45 it 10 is in perfect accord with lack of 
detections in clusters poor in Blue Stragglers where their absence can be simply explained by the 
Poisson fluctuations. 

The high frequency of contact binaries among Blue Stragglers is visible in the color index 
and period distributions shown in Figure |8[ These distributions are obviously far from being 
rigorous in the statistical sense, yet they do show interesting trends. When compared with the 
Disk systems (R98), the GC contact binaries occupy only the blue ends of the {B — V)o and 
(V — I)o distributions. Note, however, that the blue end points are the approximately same for the 
GC and Disk sample distributions indicating that the Disk sample may contain an admixture of 
low-metallicity objects similar to those in globular clusters. While the red systems are mostly likely 
under-represented in the GC sample because of the selection biases against faint, red systems, we 
see a definite lack of long-period systems which are intrinsically the brightest and should be easily 
detected. Of course temporal-window biases for the GC sample may have contributed here (in 
the sense that monitoring programs were by necessity short), especially when compared with the 
excellent data for the 5 kpc OGLE sample (R98) which defines the long-period part of the period 
distribution in Figure ^. 



11. Conclusions 

Although at least 1/3 among 86 systems presumably located in the analyzed globular clusters 
are foreground projections from the Disk which must be carefully weeded out, we can state with 
confidence that contact binaries in globular clusters are definitely different than the very common 
Disk population systems. The main feature are their short orbital periods resulting from small 
dimensions of components. This is seen not only in the distribution of the orbital periods, but also 
in low luminosities. Thus, the downward shift of the contact binary subdwarf sequence below that 
for the Disk systems is primarily due to the reduced-dimension effect, not to the blue shift caused 
by the reduced blanketing. The long-period systems are intrinsically more luminous and easier to 
discover so that their absence is highly significant. 

While the metallicity effects are clearly seen in the properties of the Population II contact 
binaries, to our surprise, they are not visible in the My"' = My'^log P, color) calibrations. More 
exactly, they do not manifest themselves in the calibration based on the (V — I) color index; the 
data utilizing the {B — V) color index are too poor to be sure that the calibration does not need a 
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metallicity term. For simplicity, we assumed that both caUbrations (which are used only to select 
the systems, not to analyze them) do not require any [Fe/H] dependent terms. It is recommended 
that the VI bandpass filter combination be used in the future because the BV data may indeed 
show some metallicity dependence, but - more importantly - are more susceptible to photometric 
errors for red stars and to reddening correction uncertainties. 

Very little can be said about contact binaries located on the Main Sequence, below the 
Turn-Off Point (TOP). Severe discovery selection effects are suspected from the peculiar 
distribution of amplitudes with missing small amplitudes. In contrast, 20 Blue Straggler contact 
binaries known at this time give a reasonable estimate of their frequency of occurrence of one such 
system (counted as one object) per 45 ± 10 Blue Stragglers. This reciprocal apparent frequency 
is about 2-3 times higher than for the Disk systems among the normal F-K stars. It is entirely 
possible that the same mechanism which produces a continuous sequence of contact binaries across 
the TOP in old open clusters simply had more time to produce more contact systems in globular 
clusters. 
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Captions to figures: 

Fig. 1. — The absolute-magnitude deviations AMy are shown here for both data-sets {BV set in the 
upper panels, VI set in the lower panels), versus the respective colour indices and logP (in days). 
Note the large fraction of foreground systems with negative values of AMy. Most systems with 
blue colors are apparently genuine members and also Blue Stragglers. The approximate locations 
of the Turn-OfF Point for the "metal-rich" GC's are shown by the thin broken vertical lines (see 
Figure 5). The Class-1 members are those with AMy < 0.5 mag.; the Class-2 members are defined 
by larger deviations, but smaller than one magnitude. 

Fig. 2. — Histograms of the deviations AMy with the values of My*' computed without any 
metallicity corrections (continuous lines) and with the metallicity corrections of —0.3 x [Fe/H] for 
the BV set and —0.12 x [Fe/H] for the VI set added to Equations || and ^ (broken lines). See the 
next figure for the details within a narrow range of the deviations. 

Fig. 3. — The band within —2 < AMy < +1 mag. is shown versus the metallicity parameter 
[Fe/H]. The clusters are identified by NGC numbers along the upper edges of the panels. The 
slanting lines give the corrections oc —0.3[Fe/H] and oc —0.12[Fe/H]. The distribution of data 
points around the zero-deviation line with and without the metallicity corrections is the main 
argument for abandoning their use in the calibrations. The data for NGC 4372 are shown by open 
circles with arrows showing the amount of shift for all systems in the cluster under an alternative 
set of assumptions on (m — M)y and Eb-v (see the text). 

Fig. 4. — Radius changes expected from metallicity variations following the study of Castellani et 
al. 1999. The parameter C = 1^ / describing the correlation of the helium abundance change 
with the heavy- metal abundance change is currently very poorly known. 

Fig. 5. — Color- magnitude diagrams for contact binary systems in globular clusters. Filled circles 
mark Class-1 member candidates which are most probably the real members of the clusters. Open 
circles mark Class-2 member candidates which are probably not members. Metallicities of individual 
systems are coded by the size of the symbol with large symbols for metal-poor systems with 
[Fe/H] < —1.5. The expected metallicity effects are shown by arrows in the left corners of each 
panel for [Fe/H] = —1 and —2. The vertical broken lines give the approximate limits for occurrence 
of Blue Straggler systems at {B — V)o = 0.4 and {V — I)o = 0.55. Small crosses show the data 
for the Disk Population contact systems in open clusters (R98); the {B — V) data have been taken 
directly from this paper whereas the {V — I) data are from the original papers cited there. 

Fig. 6. — Period-color diagrams for contact binary systems in globular clusters. The symbols used 
here are identical as in the previous figure. The Short Period Blue Envelopes (SPBE) for the Disk 
stars are shown as curved lines, as described in the text. The arrows giving the expected variations 
in the period and color due to metallicity changes for T^f / = 5000 K have the origin slightly shifted 
for better visibility from the nominal values of (B — V)o = 0.88 and (V — I)o = 0.93. To give a 
feel on how absolute magnitudes relate to the plotted quantities, the slanting dotted lines across 
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the figures give the values of My"' (not My^ as in the previous figure as these are individual to 
systems). 

Fig. 7. — Amplitudes of light variations for both color-index data-sets are shown in relation to de- 
reddened color indices. The vertical broken lines delineate the regions of Blue Stragglers. The data 
with crosses superimposed on the circles (closed for Class- 1 and open for Class-2 members) mark 
the EB-type systems. Since their primary eclipses are deeper than for normal contact systems, they 
have - on the average - larger amplitudes than typical contact system. 

Fig. 8. — The color-index distributions for Class-1 members (continuous line histograms) with 
added Class-2 members (broken line histograms) are compared with Disk systems from the Baade's 
Window OGLE sample (R98) (shaded area). Note that while the OGLE data come from a volume- 
limited samples (3 kpc, with modifications at the long-period end from the 5 kpc sample), the 
GC sample most probably contains many detection-selection biases. The OGLE sample data are 
shown here with arbitrary normalization, only to indicate the shapes of the distributions. The 
upper panels give the BV set data, the lower panels give the VI set data. 
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Table 1. Parameters of the globular clusters with contact binaries. 



Name 


Name 


I (deg) 


b (deg) 


Roc (kpc) 


Eb-v 


(m - M)v 


[Fe/H] 


c 






Nbs 


NGC 104 


47 Tuc 


305.90 


-44.89 


7.4 


0.04 


13.37 


-0.76 


2.03 




8/7/0/2 


250 


NGC 288 




152.28 


-89.38 


11.6 


0.03 


14.69 


-1.24 


0.96 


1/0/0/0 




30 


NGC 4372 




300.99 


-9.88 


7.1 


0.39 


15.01 


-2.09 


1.30 


8/0/4/1 


8/0/2/1 


30 


Rup 106 




300.89 


11.67 


18.5 


0.20 


17.25 


-1.67 


0.70 


2/0/0/0 




25 


NGC 5139 


uj Can 


309.10 


14.97 


6.4 


0.12 


13.97 


-1.62 


1.24 




18/13/2/11 


300 


NGC 5272 


M3 


42.21 


78.71 


12.2 


0.01 


15.12 


-1.57 


1.84 


1/0/0/0 




20 


NGC 5466 




42.15 


73.59 


17.2 


0.00 


16.15 


-2.22 


1.32 


2/2/0/2 




30 


NGC 5904 


M5 


3.86 


46.80 


6.2 


0.03 


14.46 


-1.29 


1.83 


0/0/0/0 


0/0/0/0 


50 


NGC 6121 


M4 


350.97 


15.97 


5.9 


0.36 


12.83 


-1.20 


1.59 


8/7/1/2 




30 


NGC 6362 




325.55 


-17.57 


5.3 


0.08 


14.79 


-0.95 


1.10 


3/0/0/0 


1/1/0/0 


25 


NGC 6397 




338.17 


-11.96 


6.0 


0.18 


12.36 


-1.95 


2.50 


2/1/1/1 




20 


NGC 6441 




353.53 


-5.01 


3.5 


0.44 


16.62 


-0.53 


1.85 




3/0/0/0 


30 


NGC 6752 




336.50 


-25.63 


5.2 


0.04 


13.13 


-1.56 


2.50 


7/0/3/1 




10 


NGC 6838 


M71 


56.74 


-4.56 


6.7 


0.25 


13.75 


-0.73 


1.15 




4/4/0/0 


50 



Notes to Table 1. 

All data in all columns with the exception of the last three are from the June 22, 1999 version of the compilation of Harris, which is avai 
http://physun.physics.mcmaster.ca/Globular.html. The columns are explained in the text. The columns labeled nP^ and nX' give the nu 



contact systems observed in both color-index sets in the format: Udet / nci / nc2 / nss ■ ndet is the total number of systems in the direction of th 
while nci and nc2 are the numbers of systems of both membership classes, as defined in the text; ubs is the number of Blue Straggler 
combined for both classes. Nbs is the approximate number of Blue Stragglers in the cluster which were monitored for variability. 

References for Table 1. 

(1) Kaluzny et al. 1998b; (2) Kaluzny et al. 1997a; (3) Kaluzny & Krzeminski 1993; (4) Kaluzny et al. 1995; (5) Kaluzny et al. 1996; (6) Ks 
al. 1997b; (7) Kaluzny et al. 1998a; (8) Mateo et al. 1990; (9) Yan & Reid 1996; (10) Drissen & Shara 1998; (11) Olech et al. 1999; (12) Kalu2 
1999; (13) Kaluzny et al. 1997c; (14) Mazur et al. 1999; (15) Rubenstein & Bailyn 1996; (16) Kaluzny 1997; (17) Layden et al. 1999; (18) T] 
et al. 1999; (19) Yan & Mateo 1994 



